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Synopsis 

Fracture mechanics is the study of the failure of a body that contains a flaw. In the 
energy balance approach to fracture mechanics, contributions from the external work and 
elastic strain energy are accounted for but rarely are corrections for the kinetic energy 
given. Under slip-stick conditions, part of the external work is expended as kinetic 
energy. The magnitude of this kinetic energy depends upon the shape of the crack. A 
specimen with a blunt crack will fail at a high load and the crack will catastrophically 
travel through the material until the kinetic energy is dissipated. Material with a sharp 
crack will fail at a lower load but will still be catastrophic in nature. 

A kinetic term is incorporated into the energy balance approach. This term accounts for 
the velocity of the crack after failure and how far the crack travels before arresting. This 
correction makes the shape of the initiation crack irrelevant. When applied to data 
generated by tapered double cantilever beam specimens under slip-stick conditions, the 
scatter in the measured critical strain energy release rate is significantly reduced. 


Introduction 

Fracture mechanics is the study of the failure of a body that contains a flaw. Griffith 
proposed that failure occurs when sufficient energy is released by growth of this flaw to 
create the new crack surface 1 . The released energy is assumed to come from the elastic 
or potential energy of the loading system and no other corrections are used. This 
approach works well for test conditions where the crack propagation is quasi-static and 
any additional crack growth results from an increase in applied load. This analysis 
breaks down when the crack propagates in a slip-stick behavior. 


This type of crack growth is characterized by a saw-tooth shaped load profile for a 
tapered double cantilever beam (TDCB) specimen as shown in Figure 1. An initiation 
load is associated with the profile as the specimen is loaded until the crack 
instantaneously propagates. As the velocity of the crack decreases and stops, an arrest 
load is reached. This type of behavior can be repeated several times as the crack moves 
down the specimen length. The load required for crack initiation is greater than that ^ 
required for stable crack growth and part of the energy is dissipated as kinetic energy . 

Under slip-stick conditions, three features complicate the fracture behavior: inertia forces, 
rate dependent material behavior, and reflected stress waves. In the case of inertia 
effects, the crack grows abruptly because a portion of the work that is applied to the 
specimen is converted to kinetic energy. Rate dependent material behavior is important 
when considering polymeric materials, because they are typically viscoelastic and very 
sensitive to temperature and deformation rates. Reflected stress waves are of concern 
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because they propagate throughout the material and reflect off the free surfaces. These 
reflected waves influence the local crack tip stress and strain field that affect the fracture 
behavior 3 . 


Maugis describes the slip-stick behavior in terms of viscoelastic losses or internal friction 
at the crack tip 4 . His model describes the critical energy release rate, G , versus the crack 
velocity, v , as shown in Figure 2. This figure is the superposition of a curve for brittle 
fracture with dynamic effects, and a curve for viscoelastic losses. He suggests that as 
long as G <G C (Figure 2) the crack propagates at constant velocity at constant G . At G c 
(Point A), the crack velocity jumps to the second curve (Point B) and the crack 
propagation is catastrophic. The velocity is too high and the crack slows down to Point C 
where the velocity is still too high, then jumps to Point D where the crack suddenly 
arrests. The velocity increases to v, and jumps to Point B. A slip-stick motion thus 

occurs. 

This approach does a good job describing the slip-stick behavior shown in Figure 1 but 
does not account for the variability in velocities observed during catastrophic failure. In 
certain experimental conditions, the magnitude of the initiation load can vary 
significantly with concomitant variation in the crack velocities. Several studies have 
examined the relationship between the crack driving force and the crack velocity 5 ' 8 . The 
magnitude of the crack driving force and hence the crack velocity strongly depends upon 
how the sample initiation crack is started. 

Initiating a proper crack can be very difficult. Marshall et al, found that for polystyrene 
the critical stress intensity factor (evaluated at crack instability), which should be a 
constant, independent of test method, and specimen geometry, varied over a wide range . 
This was shown to be dependent upon the method of notching employed. This behavior 
was also observed by Hine et al, who were able to reduce the variability in slip-stick 
behavior by reducing the probability of a region of plastic deformation occurring in front 
of the starter crack 10 . Dillard et al, describe the importance and difficulty in generating 
an appropriate initial flaw 11 . There is some uncertainty as to whether the initiation crack 
is ever adequate to give an accurate measure of the resistance of the material to fracture 
during slip-stick behavior 

In this work, we use the TDCB geometry to measure the critical energy release rate for an 
adhesive that exhibits slip-stick behavior. A kinetic energy term is included in the energy 
balance approach. The resulting expression accounts for the variability introduced by the 
different types of initiation cracks. 


Experimental 

Materials 


The epoxy resin used in this study was TIGA 321® manufactured by Resin Technology 
Group (RTG). The curing agent used in the TIGA 321® is a mixture of amines. The 
adhesive contains both Kevlar® fibers and milled glass fibers. 
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Preparation of TDCB Specimens 


The epoxy was vacuum mixed and then coated onto a D6AC steel TDCB adherend that 
conformed to ASTM D 3433. A schematic of the specimen is shown in Figure 3. Excess 
adhesive was applied to the beam and allowed to squeeze out as the two adherends were 
pressed together. The extruded adhesive was removed from the side of the beams prior to 
curing. The bondline thickness was set at 0.127 cm using Teflon® spacers controlled to a 
tolerance of 0.0127 cm. After coating, the specimens were placed in a jig that would 
maintain alignment and a slight compressive load was applied. The samples were cured 
at 40.5°C for 48 hours and then rapidly cooled to 22.2°C. This cure cycle was chosen 
because it duplicates the cure schedule of the adhesive as it is used on one of the products 
produced by the Thiokol Propulsion Group of Alliant Techsystems. 

Critical Strain Energy Measurements 


Testing was conducted using a screw-driven testing machine (e.g., Instron® or Satec®) 
outfitted with a temperature-control chamber. The tests were most often performed in 
displacement control at a crosshead speed of 0.0127 cm/minute. Unless otherwise 
specified, all tests were performed at 22.2°C. At this test temperature and deformation 
rate, TIGA 321® displayed slip-stick behavior. 

A variety of techniques were used to generate flaws in the adhesive bondline. Among 
these were embedding a razor in the adhesive bondline during fabrication, tapping a razor 
blade into the adhesive bondline at cold temperature, fatigue loading the specimen at cold 
temperature, or loading the beam and exposing the bondline front to liquid nitrogen. 

Most of these approaches introduced a visibly sharp crack into the bondline (with the 
embedded razor blade being a questionable exception). However upon loading the 
specimens at the specified temperature and deformation rate, a load would be reached at 
which the crack would propagate down the entire length of the beam. The specimens 
with the embedded razor bladed also reached a critical load but the crack would only 
propagate a short distance, allowing more than one loading peak before failing 
catastrophically. Because of this, the embedded razor was selected as the method for 
inducing a flaw into the bondline. 

The crack velocity was measured using crack propagation gauges manufacture by 
Micromeasurements (part number TK09-CPC03-003/DP). The gauges consist of parallel 
electrical circuits that produce a change in voltage, as the circuits are broken. These 
gauges are 3.86 cm in length, which is shorter than the full length of the adhesive 
bondline, but this will be addressed later in the results and discussion section. 
Measurements were taken from the gauge at a rate of 222,000 readings/s. 

Measurement of the Acoustic Velocity in Cured TIGA 321 

The measurement of the acoustic velocity in the adhesive material was done using a pulse 
multiple echo technique, which can simultaneously measure the acoustic wave-speed and 
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the thickness of the sample. A pulsed ultrasound wave from a 1MHz (center frequency) 
transducer was excited by a HV pulsed from a Ritec SP-801 pulse to the transducer. The 
beam was directed normally through one of the ends of the adhesive sample and echoes 
from the front surface, back surface, and the surface of another passive transducer located 
behind the sample were reflected back to the active transducer. This is shown in Figure 
4. A Ritec RDX-6 Diplexer was used to separate the transmitted pulse signal from the 
transducer signal produced by the reflected energy. A Ritec B40 Broadband Receiver 
amplified the electric signals generated by the reflections striking the transducer. The 
amplified signal was digitized by a Lecroy LC9314L digital oscilloscope and stored in a 
computer file for later reference and analysis. 

The theory behind the technique is explained as follows. By examining Figure 4, it can 
be seen that the difference between the arrival times of the reflector reflected waves with 
and without the sample in the water medium are 


2d 2d 



where t is the arrival time of the reflector reflected waves without the sample, t rs is the 

ro 

arrival time of the reflector reflected waves with the sample, d is the thickness of the 
sample, v, is the acoustic velocity in water, and v 2 is the acoustic velocity in the cured 

adhesive. 

The difference in arrival times of the front surface and back surface reflections are 


h l f 


2d 


v 


where t f is the arrival time from the sample front surface and t b is the arrival time from 
sample back surface. 

Thus, 


2d — (t ro t rs +t b t f V, 


and. 


v 2 = 
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t b~ t f 
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Results and Discussion 


Slip-Stick Fracture Behavior 

The dependence of the critical strain energy release rate for TIGA 321® upon temperature 
and deformation rate is shown in Table 1. Included in the table is the coefficient of 
variation (CV), which is the standard deviation, divided by the mean average. The 
reported values are evaluated at crack instability using traditional relations for calculating 
the strain energy release rate for a TDCB specimen as discussed in ASTM D 3433. 
Inspection of the table shows that the material exhibits slip-stick behavior at cold 
temperatures and high deformation rates as expected. Of particular concern is the 
considerable data scatter under these test conditions. These fracture energy values can be 
used to establish safety factors, and if a statistical penalty is applied to the mean average 
(e.g., 3-sigma), the resulting values will be extremely low. Another problem with basing 
the fracture energy on measurements taken at crack instability is that a portion of the 
calculated energy is kinetic energy. Also, the failure load has a strong dependence upon 
crack shape. 

If the slip-stick behavior is closely examined, certain observations can be made. For test 
conditions of 22.2°C and 0.0127 cm/minute, two load peaks are observed. A typical load 
profile is shown in Figure 5. The sample is loaded until a critical value is reached. The 
crack then propagates down the length of the beam in a catastrophic manner decreasing in 
velocity until the crack halts. The beam then continues to load under the constant applied 
displacement rate until a new critical load is reached and the crack propagates down the 
remaining length of the beam in a catastrophic manner. If the load of the first failure 
peak is plotted versus crack length obtained in the first unstable crack growth event as 
shown in Figure 6, a linear relationship is observed. The more interesting result is 
observed when the second peak is plotted versus crack length from the first catastrophic 
failure event as shown in Figure 6. There is again a linear relationship between the load 
and the crack length, but it behaves in an opposite manner to that of the first load peak. 
Generally, the higher the load is on the first peak, the lower the second peak will be. This 
phenomenon can be repeated by unloading the sample after reaching the first peak, 
staging the specimen for 3 days under ambient conditions, and then reloading. Even after 
the staging period, the second peak follows the same behavior indicating that this 
behavior is associated with the adhesive and not the loading setup. 

These types of trends have been previously observed in the literature. Hine et al, 
observed that the crack velocity was proportional to the crack length and also to the 
initiation load 10 . This is consistent with our observations. For dissipative processes that 
involve an increase in the energy release rate associated with crack growth, the material 
is said to exhibit R-curve behavior. The lack of an R-curve indicates that fracture is only 
associated with the intrinsic toughness of the material and no other toughening 
mechanism is involved 12 . Du et al, observed that the intrinsic toughness for their 
modified epoxy decreases with increasing crack velocity 13 . Generally, there is an inverse 
relationship between the intrinsic toughness and crack velocity 3 . 
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The implications of this inverse relationship are a bit disturbing. As discussed 
previously, the load at failure depends upon how the specimen is precracked, and the 
general consensus is that a sharp initial crack is ideal. Cracks generated by slip-stick 
behavior are extremely sharp and yet subsequent loading gives considerable scatter in the 
failure loads. Because of this, there is always an associated uncertainty as to whether the 
crack is generated under the proper conditions to give an accurate measure of the critical 
strain energy release rate. 

Kinetic Energy Contribution to the Critical Strain Energy Release Rate 

The crack shape will affect the crack velocity during catastrophic failure. A sharp crack 
will fail at lower loads with the result of a slower crack velocity. The opposite trend in 
velocity will be followed for blunt cracks. If the behavior of the velocity can be 
normalized to some standard, the effect of the crack shape becomes irrelevant. This 
correction is found by calculating the kinetic energy expended during catastrophic failure. 
In the energy balance approach to fracture mechanics, quasi-static crack propagation is 
described as the conversion of the work done by the external force and the available 
elastic energy stored in the bulk of the specimen into surface free energy. If the crack is 
not quasistatic but propagates unstably in a slip stick faction, a kinetic energy term is 
required in the energy balance. For work on a system, the total energy balance for 
external work, W , is given by 


W =U +T + D 

where U is the elastic strain energy, T is the kinetic energy, and D is the dissipative 
energy associated with fracture. 

For a body with a given cross-sectional area, A , the condition for crack growth is 

dD _d(W -U -T) (1) 

dA ~ dA 

The term on the left hand of the expression is called the critical strain energy release rate, 
G c , and is a measure of the energy required to propagate a crack over some unit area. It 
is assumed that the body is of uniform thickness, b , in which case equation (1) becomes 

1 d(W-U-T) 
c b da 


Yaniv et al., derived the critical strain energy release rate for a specially configured 
TDCB 14 . A similar derivation will be presented here for a differently configured TDCB. 
The kinetic energy for a beam is given by the Bemoulli-Euler beam theorem as 
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where A(x) is the cross-sectional area, p is the material density, a is the current length 
of the crack, y is the beam deflection, and x is the location along the beam length. 


The transverse velocity can be written by the chain rule as 

dy _ dy da dy dx _ dy . 

dt da dt dx dt da 

Mostovoy et al., gives the displacement for a single cantilever beam as 15 


y = 


12 Pr 


Eb { h 


[—,dx 

J u 3 


Skr)}!* 

Eb J 0 h 


(3) 


(4) 


where P is the applied load, E is Young's modulus of the beam, h is the height of the 
beam, and v is Poisson's ratio of the beam. 

Differentiating the displacement with respect to crack length and setting v = 1/3 gives 


dy 

da 


4 P 
Eb 


3a 2 J_ 
h? h 


AP 

- — m 

Eb 


(5) 


For a TDCB, the term is brackets is assigned to a constant, m , and is given the value 
35.43 cm' 1 . In the actual fabrication of the TDCB, the adhesive bondline doesn't start at 
the loading point but rather where the tapering of the beam begins as shown in Figure 7. 
For the bondline area, the height as a function of x is given as 

3x 2 1 

1- — = m 

/i 3 h 

The real root of this expression in terms of h can be solved. However, the solution is 
clumsy to work with numerically but can easily be Fit with a power law in the form 

h = 0.445x 0663 (6) 

A comparison of the fit between this power law and the numerical solution to the roots of 
m is shown in Figure 8. It can be seen that there is an excellent agreement between the 
two curves. 


For our experimental conditions, the load, P , given in equation (5) will not be constant 
but will vary from an initiation value, P , , to an arrest value, P A . It is assumed that the 
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load profile will be linear as the crack propagates. Thus, the load as a function of x is 
given the form 


F(*)=p,-(r,-rJ± (V 

The crack velocity, a , is also assumed to be a linear function ranging from an initial 
velocity, d„ , at crack initiation and decreasing to zero at crack arrest. The velocity as a 

function of x is assumed to be in the form 

i( x ) = ikzA ( 8 ) 

a 

Substituting equations (6)-(8) into equation (2), integrating with respect to x , and then 
differentiating with respect to a yields 

= 0.331 ^^-[P 2 +0.831 PP A + 0.369 P A 2 ]a 0663 (9) 

b da E 2 b 2 ' ' 

In equation (1), the contribution to the critical strain energy release rate from the work 
and elastic strain energy is given by 2 

i a(w'-c/) _ p i ac (10) 

b da 2b da 

where C is the compliance of the beam. Mostovy et al., gives for a TDCB with a 
constant m 15 


1 d(W-U) P 2 8 
- = m 

b da 2b Eb 


The critical energy release rate for slip-stick behavior in the TDCB can be found by 
subtracting equation (9) from equation (11). 


A plot of the — — — — ^ versus crack velocity is shown in Figure 9. In 

da / da 

generating this plot, P, — 4003 N, P t = 1334 N, b = 2.54 cm, E = 200 GPa, a = 7.62 
cm, p = 7.84 g/cm 3 . Inspection of the figure shows that the contribution from the kinetic 
energy remains insignificant until relatively large velocities are reached (e.g., 101.6 
m/sec). Once these velocities are reached, the kinetic contribution becomes more 
significant. 


Equation (9) requires the velocity and crack length to be known. Under certain test 
conditions for TIGA 321® bonded to a TDCB specimen, the crack will propagate 
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catastrophically down the length of the beam. These are considered no tests. Only test 
conditions that generated two load peaks are used. 

The crack velocity was measured using crack propagation gauges. A typical fit of the 
change in voltage with time is shown in Figure 10. This profile was consistent for all 
specimens tested. The crack velocity is linear until about half the gauge is traveled then 
the crack starts to accelerate. This phenomenon is unrelated to the kinetic energy 
generated during the failure event. Instead it is believed that this is caused by the crack 
meeting reflected waves. These waves are generated during the failure event and travel 
down the length of the beam until a free surface is reached. The reflected waves return to 
the point of initiation and collide with the advancing crack tip. 

This assumption seems reasonable if the speed of sound for TIGA 321 is compared to 
the linear component of the crack velocity. Using the acoustic emission technique 
described in the experimental section, the speed of sound for TIGA 321® was measured 
at 2794 m/s. The time it takes for the generated wave at the measured velocity to travel 
to the end and back before it encounters the advancing crack was calculated. The 
difference between the calculated and experimental time was on average 4%. This 
appears to be a reasonable explanation for the observed behavior but it is realized that the 
complete picture of the underlying mechanisms will be more complex. 

The data from a set of TDCB specimens tested at 22.2°C and 0.0127 cm/minute is given 
in Table 2. Included in the table is the initiation load, the arrest load, initial crack 
velocity, and crack length. There is considerable scatter in the initiation load, which 
indicates the variation in crack shape at initiation. The arrest load value has little scatter 
and is generally accepted to be a material property. It has been observed previously that 
the crack will arrest when the energy release rate falls approximately to that of the epoxy 
matrix 16 . However when reflected stress waves are significant, the dynamic arrest 
toughness of the material will be less than the true material resistance 17 . Because of this, 
care must be taken in assigning a material property to the arrest load. The reported initial 
velocity and crack length are generally a linear function of initiation load. 

The critical strain energy release rates obtained with and without removing the kinetic 
energy term are shown in Figure 11. The energy release rate with the kinetic energy is 
given by equation (11). The kinetic energy is corrected for by subtracting equation (9) 
from equation (11). The following values were used: b — 2.54 cm, E = 200 GPa, 
p = 7.84 cm/in 3 . The figure shows that accounting for the kinetic energy reduces the 
scatter in the data and gives a more accurate result for the material property. The scatter 
changes from a coefficient of variation of 22.8% for the uncorrected calculation to 9.6% 
for the corrected one. While the reduction in scatter is substantial if does not eliminate it. 
It is likely that a significant portion of the remaining scatter is caused by the assumption 
of a constant crack velocity. Figure 10 shows the crack beginning to accelerate after 
crack initiation. As mentioned previously, this non-linearity was mostly likely cause by 
reflected waves. Because of this, the crack will travel a different distance that it would if 
there were no reflected waves. Unfortunately, the reflected waves are an inherent part of 
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the TDCB specimen. Even with this limitation, accounting for the kinetic energy gives 
significant improvements in accuracy. 


Summary 


The general conclusion drawn from this study is that a significant portion of the scatter in 
the strain energy measured under slip-stick conditions is a result of kinetic energy. The 
magnitude of this kinetic energy contribution depends upon the shape of the crack. A 
blunt crack will cause the TDCB specimen to fail at high loads with a significant portion 
of that energy dissipated as unstable crack propagation. Sharper cracks will also cause 
unstable failure, but the velocity of the crack will be smaller in magnitude than for the 
blunt crack. 

A correction to the energy balance approach that accounted for contributions from the 
kinetic energy was presented. Using this correction, the scatter in the critical stain energy 
release rate was significantly reduced. Also, the shape of the starter crack was irrelevant 
when the correction factor was used. 
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Table 1 - Critical strain energy results for TIGA 321®. Results shown are the mean 


0 — - — 

Deformation Rate 
(cm/minute) 

Test Temperature (°C) 

GIc (J/m 2 ) 
Average (CV) 

Failure Type 

1.27 


1273 (18.5) 

Slip-Stick 

0.127 

- 6.7 

1394 (45.1) 

Slip-Stick 

0.0127 

1931 (10.1) 

Slip-Stick 

0.00127 


2326(16.4) 

Slip-Stick 

1.27 


1779 (29.8) 

Slip-Stick 

0.127 

4.4 

2304 (44.4) 1 

Slip-Stick 

0.0127 

2646 (37.4) 

Slip-Stick 

0.00127 


3366 (17.4) 

Slip-Stick 

1.27 


4789 (9.3) 

Slip-Stick 

0.127 

22.2 

5091 (9.4) 

Slip-Stick 

0.0127 

5626 (7.2) 

Stable 

0.00127 


5225 (8.9) 

Stable 

1.27 


6171 (2.5) 

Stable 

0.127 

32.2 

5964 (3.0) 

Stable 

0.0127 

5861 (4.9) 

Stable 

0.00127 


5792 (3.5) 

Stable 

1.27 


6486 (5.1) 

Stable 

0.127 

46.1 

7033 (2.8) 

Stable 

0.0127 

6311(6.2) 

Stable 

0.00127 


4505 (5.2) 

Stable 
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Table 2 - Data from r 

[DCB Specimen Testing 

Initiation Load 

Arrest Load 

Initial Crack Velocity 

Crack Length 

(N) 

(N) 

(m/s) 

cm 

2740 

1477 

116 

4.3 

3465 

1446 

183 

6.6 

3220 1 

1481 1 

163 

5.3 

3727 

1419 

204 

8.1 

4039 1 

1361 

229 

9^6 

2758 ‘ 1 

1521 

120 

3.8 

3683 

1477 

200 

8.9 

3505 

14777 

185 

9.1 

6269 

1468 

164 

7.4 

3065 

1530 

146 

6.1 

3487 

1508 

183 

7.9 

3692 

1446 

201 

9.1 

3432 

1334 

222 

1 11.2 

3821 

1419 

212 

9.9 

2896 

1570 

132 

5.8 

2740 

1477 

115 

4.3 
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Figure 1 : Slip-stick load profile for a TDCB specimen. 
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Figure 2: Theory proposed by Maugis to explain slip-stick behavior. 
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Figure 3: Bonding setup for TDCB specimen. 
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Figure 5: Typical load profile for TIGA 321® loaded at 0.0127 cm/min at 22.2°C. 
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Figure 6: First and second load peaks versus crack length for TIGA 321® loaded at 

0.0127 cm/min at 22.2°F. 
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Figure 7: Loading schematic forTDCB specimen. 
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Figure 10: Plot of voltage versus time for crack propagation gauges 







